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ABSTRACT 
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, Context. The VIMOS VLT Deep Survey (VVDS) was undertaken to map tiie evolution of galaxies, large scale structures, and active galaxy 
nuclei from the redshift spectroscopic measurements of ~10^ objects down to an apparent magnitude Iab = 24, in combination with a multi- 
wavelength acquisition for radio, infrared, optical, ultraviolet, and X-rays data. 

Aims. We present the evolution of the comoving star formation rate (SFR) density in the redshift range < z < 5 using the first epoch data 
O release of the VVDS, that is 1 1564 spectra over 2200 arcmin^ in two fields of view, the VVDS-0226-04 and the VVDS-CDFS-0332-27, and 
^ , the cosmological parameters {Cl^, Qa, h) = (0.3, 0.7, 0.7). 

C/2 ■ Metiiods. We study the multi-wavelength non dust-corrected luminosity densities at < z < 2 from the rest-frame far ultraviolet to the optical 

■ passbands, and the rest-frame 1500 A luminosity functions and densities at 2.7 < z < 5. 
> Results. They evolve from z = 1.2 to z = 0.05 according to (1 + zY with x = 2.05, 1.94, 1.92, 1.14, 0.73, 0.42, and 0.30 in the FUV-1500, NUV- 
k> ' U-3600, B-4400, V-5500, R-6500, and 1-7900 passbands, respectively. From z = 1.2 to z = 0.2 the S-band density for the irregular-like 

rS ' galaxies decreases markedly by a factor 3.5 while it increases by a factor 1.7 for the elliptical-like galaxies. We identify several SFR periods; 
from z = 5 to 3.4 the FUV-band density increases by at most 0.5 dex, from z = 3.4 to 1.2 it decreases by 0.08 dex, from z = 1.2 to z = 0.05 it 
declines steadily by 0.6 dex. For the most luminous Mab(1500 A) < -21 galaxies the FUV-band density drops by 2 dex from z = 3.9 to z = 1.2, 
and for the intermediate -21 < Mas(1500 A) < -20 galaxies it drops by 2 dex from z = 0.2 to z = 0. Comparing with dust corrected surveys, 
at 0.4 S z S 2 the FUV seems obscured by a constant factor of ~ 1.8-2 mag, while at z < 0.5 it seems progressively less obscured by up to 
~ 0.9 - 1 mag when the dust-deficient early-type population is increasingly dominating the 5-band density. 

Conclusions. The VVDS results agree with a downsizing picture where the most luminous sources cease to efficiently produce new stars 12 
Gyrs ago (at z - 4), while intermediate luminosity sources keep producing stars until 2.5 Gyrs ago (at z - 0.2). A modest contribution of dry 
mergers and morphologies evolving towards early-type galaxies might contribute to increase the number density of the bright early types at 
z < 1.5. Our observed SFR density is not in agreement with a continuous smooth decrease since z ~ 4. 
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1. Introduction 
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* Based on observations collected at the European Southern 
Observatory Very Large Telescope, Paranal, Chile, program 070.A- 
9007(A), and on data obtained at the Canada-France-Hawaii 
Telescope, operated by the Institut National des Sciences de I'Univers 
of the Centre National de la Recherche Scientifique of France, the 
National Research Council of Canada, and the University of Hawaii. 



1.1. The global star-forivation history one decade ago 

The history of the comoving space density of the star formation 
rate (SFR) is a key study which has undergone a sp ectacular ex 



plosion of publications since the pioneer work of iMadau et al 



(199 6) in which ga laxy surveys were used for the first time. 
Mad au et al.l ( Il996h derived the global star formation as a func- 
tion of redshift combining emissivi ties from three distin ct sur- 



veys: the Ha-selected UCM survey dGallego et aUll995l) at z < 



0.05, the largest spectroscopic sample at that time, the /-band 
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selected Canada- France Redshift Survey at 0.2 < z < 1 (CFRS; 



Lilly et al.L 11996b . and their own galaxy sample at 2 < z < 4.5 



using the Lyman co ntinuum break colour-selection technique 
( Steidel et al. , 19961) applied on the deep optical imaging, the 
Hubble Deep Field (HDF) survey. Nevertheless high-z sources 
were lower limits, and incompleteness corrections and uniform 
limiting magnitude were required. 

The SFR density is usually derived from the mean lutni- 

poo 

nosity density, or emissivity, that is X, - (f>(L)LdL. At first 
sight the SFR density appears a simple and useful tool to trace 
back the evolution of star formation and to hnk it with the evo- 
lution of stellar mass, but a large spread between different mea- 
surements have led to controversy. Uncertainties in conversion 
factors from luminosity to SFR, and in the amount of dust ob- 
scured SFR, coupled with the different selection criteria of each 
survey and with the uncertainty about the shape of the lumino- 
sity function, meant that the SFR history of the Universe was 
poorly determined, and thus hotly debated. In the mid-90's the 
main questions that arose were: is there a peak of the cosmic 
star-formation history at 1.3 < z < 2.7? By how much does 
interstellar dust attenuate the ionizing flux? Is the evolution so 
rapid below z - 1 ? Is the high-redshift dropout population rep- 
resentative? Does the red galaxy population evolve passively or 
not? 



1.2. A wealth of multi-wavelength studies 

Measuring the cosmic star-formation history has advanced 
rapidly in the last decade meanwhile one had to translate results 
using the standard cold dark matter (CDM) into the ACDM 
cosmology. It adds uncertainties since the cosmological cons- 
tant introduces a redshift dependence to the luminosity function 
unless one can take back all the data needed to compute the 
luminosity function. Still the infeiTed larger distances and vo- 
lumes result in decreasing luminosities and densities, and thus 
in a shallower SFR evolution. 

On the one hand, a wealth of multi-wavelength emissivi- 
ties (far- and near-ultraviolet, far- and near-infrared, radio. Ha, 
etc.) has been obtained, all of them being more or less di- 
rectly proportional to the ionizing ultraviolet stellar spectra at 
A < 912 A mainly produced by massive, young, short-lived 
stars (OB stars, t < few x 10^ yr). On the other hand, super- 
novae events and their by-products (neutrinos and gamma rays) 
are proportional to the star-death rate and are also used to probe 
the SFR. Nevertheless these results depend most strongly on 
stars more massive than 3 solar mass, and require extrapolation 
of the initial mass function (IMF) to lower masses to obtain the 
SFR for all masses. Futhermore, the ultimate desired value is 
the amount of the interstellar gas mass transformed into stars, 
and the correlated gas mass return into the interstellar medium 
as stars deplete their initial fuel sources; it requires assumptions 
and models for stellar atmospheres and stellar evolution tracks. 

An example of an extensive compilation drawn from the lit- 
erature o f SFR d ensit y measurements at < z < 6 w as done by 
HopkinsI (120041) . and iHopkins and BeacomI (120061) . and it led 
to a cosmic SFR history constrained to within factors of about 
~ 3. There is now growing evidence that the evolution has no 



peak at z ~ 1.5, however it is still unclear whether at z > 3 the 
evolution flattens, or declines or continues to increase. Those 
remaining uncertainties show the necessity to further investi- 
gate galaxy redshift surveys to constrain galaxy assemblies. In 
particular, one needs to explore multi-wavelength datasets over 
the same field of view, and to obtain homogeneous datasets 
over a large redshift range. 

A detailed picture of the SFR history is emerging with 
the advent of large, deep redshift surveys, coupled with multi- 
wavelength ground and space observations over the same sky 
area. These new multi-wavelength surveys are providing a bet- 
ter understanding of the nature and the evolution of the galaxy 
population (i.e.. 



e-g- 



Bell et al.L 12004 and references within). 



Indeed on one hand, one can probe the same galaxy popula- 
tion at different wavelengths which gives insights about stellar 
masses (near-infrared) and star formation (far-ultraviolet, far- 
infrared), possibly coupled with spectroscopic indexes. On the 
other hand, one can measure reliable comoving volumes with 
accurate redshifts and intrinsic luminosities of objects which 
both are key measurements to estimate luminosity functions. 
Futhermore observing in various windows of the electromagne- 
tic spectrum enable to select the largest galaxy sample which 
includes sources with specific energy distributions, and de- 
tectable at only some wavelengths. The accuracy in the mea- 
surement of the shape of the luminosity function, coupled with 
the stellar mass-to-light ratio of galaxies, is crucial to estimate 
the amount of baryons in stars at a given epoch of the Universe. 

1.3. The detailed study of the VVDS 

The Visible Multi-Object Spectrograph (VIMOS) installed 
on the European Southern Observatory (ESO) Very Large 
Telescope (VLT) was built to produce systematic large red- 
shift surveys thanks to its high multiplex capabilities, e.g. -550 
R ^ 230 s pectra of sources ob served simultaneously over 
218 arcmin^ dLe Fevre et al.Ll2003h . We conducted the VIMOS 
VLT Deep Survey (VVDS), a major multi-wavelength spectro- 
scopic survey, to investigate the evolution of galaxies. Active 
Galaxy Nuclei (AGN), and large-scale structure^ 

Here we detail our work in using the /-selected VVDS first 
epoch data described in Le Fevre et al. (2005a) to study the 
evolution of the luminosity density within the redshift range 
< z < 5. It is part of a series of papers which analyse different 
aspects r elated to th e lumin osity function evolution with this 
data set. Ilbert et al. i (l200l de scribe the global o ptical lumi- 
nosity functions over < z < 2. IZucca et al. I(l2006l) explor e the 
color-ty pe luminosity functions over < z < 1.5. Ilbert et al.l 
(l2006a ) investigate the contribution of different morph ologi- 
cal types to the luminosity functions. 'il bert et"aL I (l2006b l) anal- 
yse the lumi nosity functions in different environments over 
< z < 1.5. |PaltanieTaD(l2007h ar ialyse in detail the 17 00 A 
luminosity function at 3 < z < 4. lArnouts et al.l (l2005b des- 
cribe the 1500 A luminosity functions over < z < 1.3 using 
GALEX-VVDS data. The infrared luminosity functions with 
K-band data and with SWIRE-VVDS data, and the luminosity 
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functions for different spectroscopic-based classes are in prepa- 
ration. 

This paper is organized as follows. In Section 2 we present 
the data. In Section 3 we detail the methods used to estimate the 
comoving luminosity densities and associated uncertainties. In 
Section 4 we present the multi-wavelength global luminosity 
densities at < z < 2 derived in the UBVRI and in the near- 
UV, far-UV passbands. In Section 5 we compare our results to 
other surveys at z < 2. In Section 6 we investigate the lumi- 
nosity densities for different galaxy types at z < 2. In Section 7 
we measure the rest-frame 1500 A luminosity functions and 
densities at 2.7 < z < 5 and we compare our results to other 
surveys. In Section 8 we detail the evolution of the global far- 
UV luminosity density all the way from z = 5 to z = 0, and 
we analyse its dependency to the luminosity. In Section 9 we 
derive the history of the star formation rate density since z = 5 
and discuss the issue of dust obscuration. Finally in Section 10 
we recap our conclusions about the evolution of the rest-frame 
luminosity densities in a well controlled and homogeneous I- 
selected population over the large redshift range < z < 5 
as observed by the VVDS. Throug hout this paper we use the 
AB flux normalization (Oy. ll974 . We adopt the set {Q-u, ^a, 
h) - (0.3, 0.7, 0.7) for the cosmological parameters. 

2. Data 

Our studied sample is taken from the /-band selected spec- 
troscopic data of the first epoch observations obtained in two 
fields of view, VVDS-0226-04 and VVDS-CDFS (i.e. VVDS- 
0332-27 in the Chand ra Deep Field-South) and described in 
Le Fevre et al.l (l2005al) and lLe Fevre et al.l (120041) . It consists of 
1 1564 spectra and it covers 2200 arcmin^ of sk y area observed 
in fiv e optical passbands [/, B, V, R, and / ( Le Fevre et al.L 
20041) . In this paper we consider the well-defined selection 
function of the spectroscopic targets selected from the VVDS 
photometric parent catalogue with apparent magnitudes in the 
range 17.50 < Iab ^ 24.0. We do not use any serendipi- 
tous sources observed randomly in the slit other than the tar- 
get. Spectroscopic observations were efliciently targeted u- 
sing the VMMPS tool developed by our team for the spectro- 



graph VIMOS-VLT/ESO (see Bottini etal.Ll2005L for details) 



We used the red grism (5500 to 9500 A) and a resolution of 
R = 227. No pre-selection has been applied in terms of colors, 
sizes, photometric redshifts, or peculiar sources. The VVDS is 
based on the sole criterion of a Iab flux limit. 

Spectroscopic observations have been automatically 
processed using the VIPGI tool that we developed 



( Scodeggio et aL ,_ 2006h an d spectroscopic determination 



is described in Le Fevre et al ](l2005a). The 1 cr accuracy of the 
redshift measurements is estimated at 0.0009 from repeated 
VVDS observations. We emphasize that we have obtained 
an excellent efficiency for determining redshifts at z < 2 
and at z > 2.7. At 2 < z < 2.7 reliable spectral features are 
difficult to detect, and thus the efficiency to measure a redshift 
in this range is very poor. Observations extending further to 
the blue or into the near-IR are required to fill in this gap 
with more redshifts. We therefore present our measurements 
in the redshift ranges [0.05 - 2] and [2.7 - 5]. In total our 



targets have been classified as 7840 galaxies, 751 stars, and 
71 quasars with a reliable spectroscopic identification at a 
confidence level higher than 81 percent (coiTesponding to 
the VVDS quality flags 2, 3, 4, and 9), 1580 spectra with an 
uncertain spectroscopic identification at a confidence level 
within [48-58] percent (corresponding to the VVDS quality 
flags 1) and 792 spectra not identified (corresponding to the 
VVDS quality flags 0) . There are 7631 (1182), 31 (47) and 
178 (271) reliable (uncertain) galaxy redshifts at < z < 2, 
2 < z < 2.7, and 2.7 < z < 5 respectively. In our study, we 
exclude the quasars, which are easily identified thanks to the 
presence of large broad spectrocopic emission lines. Under the 
term 'galaxies' we note that we include any naiTow emission 
line AGN. 

As we do not have a measured redshift for every source 
to a fixed magnitude limit in the observed field of view, we 
introduced a statistical weight, which is a function of ap- 
parent magnitude and redshift and connects for sources not 
observed (Target Sampling Rate; TSR) and for sources for 
which the spectroscopic identification failed or is uncertain 
(Spectroscopic Success Rate; SSR). This statistical weight has 
be en applied to each measured galaxy at < z < 2 as described 
in lllbertet al.l(l2005h . The SSR was estimated in two ways; u- 
sing the photometric redshifts, and using the uncertain redshifts 
with a confidence level within [48-58] percent and assuming 
that the failed identifications have the same redshift distribu- 
tion. The two SSR estimates were discrepant in the redshift bin 
1 .5 < z < 2 only, with the former being twice as large as the 
latter. Using the deeper multi-wavelength observations of the 
Canada-France-Hawaii Legacy Surve}|3 over the VVDS-0226 
field, Ilbert et al. (2006c) obtained better photometric redshifts 
than in lllbertetal.1 (12005) . And thus we have refined our SSR 
estimations. In the redshift bin 1.4 < z < 2 the SSR is changed 
by a multiplicative factor of 0.50, otherwise the SSR estima- 
tions did not change within t he redshift range < z < 1 .4 with 
respect to those presented in l llbertet al.l(l2005b . 

Absolute magnitude measurements are optimized accoun- 
ting for the full information given by the multi-band photo- 
metric data in a way which minimizes the dependency on 
the templates used to fit the observed colors (see Fig. A.l in 
lllbert et al.[ l2005h . That is, we automatically choose the ob- 
served apparent magnitude which is as close as possible to the 
rest-frame band redshifted in the observer frame, so the de- 
pendency to the template is null or the smallest possible. We 
use the tem plates generated with the galaxy evo lution model 
PEGASE.2 dFioc and Rocca-Voknerangei Il997h . Finally, we 
use a sample of galaxies which are equally visible, that is 
within a given absolute magnitude range which depends on the 
rest-frame wavelength as we describe in the next section. 

3. Measuring comoving luminosity densities 

3.1. Definition 

Comoving luminosity densities, X = Jg (/i(L)LdL, depend on 
the shape of the luminosity function, ^(L)dL. In our present 
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study, the galaxy luminosity function (LF) follows a lSchechtei 
( Il976h function characterized by a luminosity, L*, a faint-end 
slope, a, and a normalization density parameter, (/>*, and thus 
X - r{a + 2). The LF is a fundamental measurement of 
the statistical properties of the population of galaxies; it is the 
distribution of the comoving number density of galaxies as a 
function of their intrinsic luminosity at a given epoch. Despite 
its simple definition, its estimation requires careful analyses of 
the survey strategy, the selection criteria, and the completeness. 

The faint-end slope is often measured in the range -1 < 
a < -2; thus for a non-diverging density of galaxies, the LF 
must have a cut-off at faint luminosities. Such a cut-off has not 
yet been observed. This implies a high-space density of low- 
luminosity galaxies, but although these galaxies are very nu- 
merous, they contribute little to the mean luminosity density; 
for instance, sources fainter than O.IL* contribute less than 20 
percent to X. for a < -1.3 (see Fig.[T]i. As we do not observe 
the faintest galaxi es, we use the Schechter functional form for 
the STY estimate dSandage et al.L ,1979.) to suppose the behav- 
ior of the LF at low luminosities. As the three Schechter para- 
meters are highly correlated, it is necessary to build the lumi- 
nosity function over the largest possible range of luminosities. 
Indeed, a weak constraint of the slope may have a strong im- 
pact on the determination of L*, which is directly translated 
into the luminosity density estimation. Nevertheless, the latter 
is more robust than the estimation of each single parameter of 
the Schechter function alone. 

3.2. Method 

We estimated the LF parameters, a, M* (or L*), and (p*, u- 
sing the Algorithm for Luminosity Functions (ALF) developed 
within the VVDS consortium. ALF uses the non-parametric 
Vmax, SWML and and the parame tric STY l uminosity func- 



tion estimators (see Appendixes in Il bert et all 2005, and refe- 
rences within). Each estimator presents advantages and draw- 
backs, and each one is affected differently by different visibility 
limits for the various galaxy types detected in deep flux-limited 
surveys. Galaxies are not equally visible in the same absolute 
magnitude range mainly due to the spectral type depe ndency 



on the k corrections. The bias was quantified in Ilbert et al 



(120041) ; it affects the faint-end slope of the global LF which 
can be over/underestimated depending on the adopted estima- 
tor. When the differences between the estimators are larger than 
the statistical uncertainties, it indicates the presence of a signif- 
icant bias. Thus in a given redshift range our LF parameters are 
estimated with data restricted to the absolute magnitude range 
in which all galaxy types are visible. It enables us to calculate 
an unbiased LF slope. 

We apply two approaches to derive the luminosity density 
from the LF. First, we derive minimal comoving luminosity 
densities in summing the LFs over the observed luminosities; 
in this case there is no assumption made over the bright or faint 
end of the LF which may be not observable in the lowest or 
highest redshift bins. Secondly, we derive the estimated mean 
comoving luminosity densities in summing the LFs over all lu- 
minosities. As we do not observe the faintest galaxies, the lat- 



ter estimates are derived by extrapolating the LF obtained us- 
ing the STY estimator. This approach is the only way to com- 
pare data through cosmic epochs since we integrate to the same 
faintest luminosity, as long as a cut-off at faint luminosities is 
not observed. 



3.3. Uncertainties 

The LF parameters are correlated to each other, and the effect of 
this correlation is that the uncertainty in the LF integral requires 
the incorporation of the LF parameter error ellipse, in addition 
to the Poisson uncertainties typically quoted. 

Thus for the two correlated parameters, a - M*, the uncer- 
tainty of X is derived from the optimal confidence regions as 
determined using the STY errors. The error ellipse implies that 
the uncertainty is not given by the squares of the one-standard- 
deviation (cr) errors of a and M* as done in the case of two 
individual, normally-distributed parameters which give an esti- 
mated value within a 68 percent confidence interval. Actually, 
the correlation increases the errors on the other parameter. 
Indeed, the probability that a and M* simultaneously take on 
values with the one-cr likelihood contour is 39 percent only. 
We use the likelihood contour corresponding to a 68 percent 
confidence interval, that is at 2.3cr. Using the number counts, 
we derive <p* for each point of this contour. Our final uncer- 
tainties coiTespond to the two points of the eiTor contour which 
give the lowest and highest values of X. 

In this procedure none of the two parameters, a or M*, has 
been fixed. If we fix one of the two parameters, then it decreases 
the errors on the other parameter. In this case, errors are usually 
not realistic, so we use another method to give uncertainties as 
follows. We derive L* T(a + 2) using the estimated Schechter 
parameters with the low, high, and mean values chosen for the 
fixed parameter, and with the Icr error on the other free para- 
meter. Using the number counts, we derive (p* for each point of 
the single free parameter axis. In the case of a fixed parameter, 
we take the extreme uncertainties given by the highest (low- 
est) X estimation calculated with the high (low) value of the 
fixed parameter, and subtracted to the lowest (highest) estima- 
tion calculated with the mean value. When we fix either M* or 
a, we footnote it in our tables of X values. 



3.3.1 . Fixing tine Scliecliter parameter M* 

The brightest galaxies are not sampled in the first redshift 

bin (0.05 < z < 0.2) because of the VVDS bright Hmit at 

Iab = 17.5. Thus in this redshift range we fix M* ± 0.05 to its 

mean local value, as found by the SDSS at z = 0. 1 in the U, B, 

J — — ' ' ' 

V, R, and I bands jBlanton et al.. 2003). We ga ve these values 

in Section 5.1 and Table 1 of Hlbert et all (l2005l). In the f ar-UV 
band we fix M* +0. 10 to the one derived by lWyder et al.l (l2005b 
from the GALEX-2dFGRS survey at < z >= 0.055 The result- 
ing STY parameters are with Mg^ = -18.12 mag, a - -1.13, 
and 0* = 0.00732 Mpc-^ 



L. Tresse et al.: The cosmic star formation rate evolution from z = 5 to z = from the VIMOS VLT Deep Survey 



5 



e 3: 

3 s 0.5 

" 

5 100 
T.S 80 
I I 60 

5 I 40 

S 

1 I 20 



(a) L*, ip' fixed 



— I 1 1 1 1 1 1 1 1 1 1 1 1 h- 

(b) L < L" 

L < 0.4L" - ■ --^ / 

L < O IL' ... " ^''^ / 



-1 -1.5 
at (LF faint end slope) 



Fig. 1. (a). Quantitative change in log of the comoving luminos- 
ity density, X, as a function of the faint-end slope, a, of the LF 
for fixed L* and <p* values, compared to Log X for a = -1. (b) 
Fraction of X as a function of a due to sources fainter than L* 
(dotted line), 0.4L* (dashed line), and O.IL* (dot-dashed line). 

3.3.2. Fixing tine Scliecliter parameter a 

At high redshift (z > 1 .2), we need to fix a because the VVDS 
faint limit, I^b = 24, prevents us from observing luminosities 
faint enough to constrain the LF faint-end slope. One of the 
most important things done to obtain a reliable estimate of X is 
to sample the LF around the knee of the luminosity distribution. 
If we do not observe magnitudes fainter than the LF knee, (i.e. 
L < L*), then X is underestimated whatever the value given to 
a. As long as a is greater than -1 the total luminosity density 
is little dependent on the slope, as illustrated in Fig. [1] and is 
dominated by the product of (p* and L*. For slopes < -1, and 
constant (p* x L*, X increases of 32, 78, and 209 percent from 
a = -1 to -L3, -1.5, and -L7, respectively. Thus uncertain- 
ties on X can be quite large when a < -1.3. The low, high, and 
mean values for a fixed slope are described in the appropriate 
sections. 

4. Global luminosity densities at < z < 2 

4.1. The VVDS multi-wavelength data 



The S chechter parameters are those derived in lllbert et al 
( l2005b in the rest-frame U-3600, B-4400, V-5500, R-6500, and 
1-7900 passbands from z - 0.05 to z = 1 . In the redshift bins 
LO < z < 1.3 and 1.3 < z < 2.0 we were very cautious in 
fixing the LF faint-end slope to the value of the slope esti- 
mated in the redshift bin 0.8 < z < 1.0. Now, we estimate 
the slope in the redshift bin 1 . < z < 1 .2, and at z > 1.2 



we use our result from Fig. 9 in Ilbert et al.l (l2005h illustrating 
the steepening of the LF faint-end slope with increasing red- 
shift. It corresponds to Aa = 0.3 between z - 0.1 and z - 0.9, 
that is an increase of 0.0375 at each step of 0. 1 in redshift. We 
have fixed a in keeping this increment since there is no reason 
why it should suddenly stop increasing. In particular, for the 



bin 1.2 < z < 1.4 (< z >= 1.3), we take the value of a mea- 
sured in the bin 1.0 < z < 1.2 (< z >= 1.1) and we increase it 
by 2 X 0.0375 = 0.075. For the bin 1.4 < z < 2.0 (< z >= 1.5), 
we take the value of a in the bin 1.2 < z < 1.4 (< z >= 1.3) 
and we increase it by 0.075 also. In these two redshift bins, we 
assume +0.075 for the error bars of the fixed a values. 

We have derived the rest-frame FUV-1500 and NUV-2800 
luminosity densities. While the rest-frame 2800 A luminosities 
are sampled by the observed reference bands from [/-band to 
/-band from z - 0.3 to z = 2.0, the rest-frame 1500 A lu- 
minosities are sampled only at z > 1 .4. Our LFs computed at 
1500 A therefore strongly rely on the spectral eneygy distr ibu- 
tion (S EP) fit at longer wavelengths. However, in .Arnouts et al.l 
(l2005h we hav e already de r ived t he 1500 A LFs using data 
from GALEX ('Martin et al.', '2005V These data were a {NUV- 
2300+/)-selected GALEX-VVDS sample from z = 0.2 to 
z = 1.2 and they were mainly based on a one-to-one identi- 
fication between optical images and NUV images. Our present 
FUV results are SED-dependent but derived from a pure I- 
selected sample similar to what was done in the other bands. 
At 0.2 < z < 0.8 the /-selected VVDS and the (NUV + /)- 
selected GALEX-VVDS LFs are in agreement within the error 
bars. This implies that even though our FUV result is SED de- 
pendent, our SED fitting is globally correct using a fixed -1.6 
slope at < z < 2. At 0.8 < z < 1.2, the two FUV luminosity 
densities differ significantly of ~0.2 dex because of a lower (/>* 
for the (NUV + /)-selected sample. The weights we adopted in 
Arnouts et al. I (12005 ) were preliminary and slightly underesti- 
mated at 0.8 < z < 1.2. 

Values of X for each passband in various redshift bins are 
given in Tab.[Tl and they are displayed in Figs.|2] |3] and |4]for 
the B, V&R&I, and U&NUV&FUV bands, respectively. 

4.2. Evolution according the rest-frame passband 

The dependence of the emissivity of the global population on 
the rest-frame band is noticeable. The non dust-corrected lu- 
minosity densities evolve with redshift over 0.05 < z < 1 .2, as 
X oc (l-Hz)" with jc = 2.05, 1.94, 1.92, 1.14, 0.73, 0.42, and 0.30 
in FUV-1500, NUV-2800, U-3600, B-4400, V-5500, R-6500, 
and 1-7900 passbands, respectively (see Figs.|2] [3] ID and |5]l. 
There is a clear differential, wavelength-dependent evolution 
of the whole population. Indeed, the average (FUV-I) rest- 
frame color emissivity of the whole galaxy population becomes 
four times redder from z = 1.2 to nowadays. Futhermore, the 
data suggest an up-turn in the emissivity evolution at redder 
wavelengths than the /-band. Nevertheless a possible up-turn 
is likely a selection effect due to the fact that the rest-frame 
/-band emissivities from an /-selected sample could be under- 
estimated in missing the very red galaxies at z > 1 . 

The FUV is related to the formation of young, massive, 
short-lived, hot stars, while the NIR is related to long-lived, 
old stars which relate closely to the stellar mass of a galaxy. 
And thus, over the last 8.5 Gyrs (z < 1.2) there has been a 
substantial decline of the star formation rate, while the old stel- 
lar conten^_showsa_smoother change in terms of emissivity. 
In iPozzetti et al.l (l2007h . we derive the stellar masses using a 
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Table 1. Comoving multi-wavelength non dust-corrected luminosity densities at < z < 2 of the VVDS with the cosmology 
(Qm, ^a, h) = (0.3,0.7,0.7). 



< z >° 


Redshift range 






log-CjsOO A 


log-C3600 A 


log -£4400 A 

W Hz ' Mpc-3 


log-Cssoo A 


log-Cesoo A 


log-C7900 A 


Minimal 


0.14 


[0.05-0.20] 


2.5 


18.65 


18.87 


19.17 


19.62 


19.77 


19.94 


20.05 


0.30 


[0.20-0.40] 


2.0 


18.83 


19.07 


19.31 


19.72 


19.87 


19.99 


20.09 


0.51 


[0.40-0.60] 


1.3 


18.82 


19.11 


19.42 


19.72 


19.91 


20.00 


20.08 


0.69 


[0.60-0.80] 


1.1 


18.99 


19.30 


19.48 


19.84 


19.98 


20.05 


20.13 


0.90 


[0.80-1.00] 


0.9 


19.00 


19.29 


19.47 


19.82 


19.95 


20.02 


20.09 


1.09 


[1.00-1.20] 


0.7 


18.97 


19.33 


19.46 


19.79 


19.92 


19.99 


20.07 


1.29 


[1.20-1.40] 


0.6 


18.91 


19.22 


19.37 


19.68 


19.78 


19.83 


19.92 


1.55 


[1.40-2.00] 


1.3 


18.81 


19.23 


19.36 


19.64 


19.73 


19.78 


19.83 


Estimated 


0.14 
0.30 
0.51 
0.69 
0.90 
1.09 
1.29 
1.55 


[0.05-0.20]" 

[0.20-0.40]' 

[0.40-0.60]'''' 

[0.60-0.80]''' 

[0.80-1.00]''' 

[1.00- 1.20]'' 

[1.20-1.40]''''' 

[1.40-2.00]''''' 


2.5 
2.0 
1.3 
1.1 
0.9 
0.7 
0.6 
1.3 


18.81-JJ5 
18.96-j^ 
19 00— " 
19 17— 

^'^■'^ '-0.09 

19 41+0 16 
19 35—"* 
19 31—^2 


18.98-;S^ 
19.10-«^ 
19 17-"2 

^^■'^ '-0.06 
19 33+0.06 
17. J J ,,03 

19 48+0 13 
19 47—" 
19 53-32 
19 65— '"' 


1 Q 9^+0.06 

19 34+0 03 

^^■-'^-0.03 

19 47+0.03 

^'■^'-0.03 

19 57-"^ 
19 69- 21 

1 Q 7^+0.06 

1 9 74+0.09 

'^-0.07 

19 74+0.13 


19 66- "'^ 

17.UU 

19 76— °'' 

'"-0.04 
1 Q 77+0.04 
' '-0.03 

19.87-;»2 

1 Q Q9+0.03 
17.7Z,_Q02 
1 Q Q7+0.03 
-''^■'^'-0.02 

19 87+0 09 

^^■°'-0.04 
19 74+0.10 
'^-0.05 


1 Q QT+0.04 
I7.7J_Q0^ 

19 96+0.09 

1 Q Q7+0.04 
^'■"-0.04 

20 03+'"'2 

^"■"-'-0.02 

20 07+" °'' 

^"■"'-0.03 

20 13+"°'' 

^"■^-'-0.04 

20 06+" '" 
20.00-;[] 


20 07+" °'' 

^"■"'-0.07 

20 10+" '2 
20 07+" °^ 

^"■"'-0.04 

20 11+" °2 

^"■^-^-0.02 

20 l4+"°5 

20.18-;?' 
20.14-;2 


20 17+"°7 

^"■^ '-0.07 

20.27+552? 
20 16+"°' 
20 19+" °3 

iU. i7_Q Q2 

90 9^+0.06 
^"■^-'-0.04 

20 27+" °' 

^"■"^'-0.05 

20 26+" '^'' 
^"■^"-0.12 

20 27+" 2" 

^"■^'-0.17 



"Mean redshift of the galaxy redshifts in the quoted redshift range. 
'Elapsed time in Gyr in the redshift bin. 
"At 0.05 < z < 0. 20 with < z >= 0.1 4 the STY LP estimators are derived using a fixed M* in the U, B, V, R, and I bands taken from the 
SDSS (< z >~ 0.1 iBlanton et al.l j2003h ; see also Section 5.1 in Ilbert et al. ( 2005)), in the FUV-1500 band from the GALEX-2dFGRS survey 



(< z >= 0.055 IWvder et al.Ll2005h . and in the NUV-2800 from our value estimated at 0.1 < z < 0.3 (i.e. -18.61 mag) since there is no local 
value in the literature. Error bars are derived from the STY estimations at M* ± 0. 10 for the FUV and NUV bands, otherwise at M* ± 0.05. 
''At 1.2 < z < 2.0 the STY LF estimators are constrained using a fixed a and error bars are derived from the STY estimations at a ± 0.0375 as 
detailed in Section 4 except for the FUV and NUV bands (see below). 

"In the NUV the slope is fixed to the value estimated at 0.2 < z < 0.4, i.e. a = -1.32. Error bars are derived from the STY high estimations 
assuming an increase of a, of -1-0.0375 each 0. 1 redshift bin, as observed in the visible (see Section 4), and from the STY low estimations 
assuming o- = -1 at 0.8 < z < 2, and from the minimal value at 0.4 < z < 0.8 since in this redshift range the latter is higher than the STY 

estimate assuming a = -I. 

'^In th e FUV the slope is fixed to o- = -1.60 at 0.2 < z < 2 corresponding to the slope found with the GALEX-VVDS sample jAmouts et all 
l2ooa . Error bars are derived from the STY estimations assuming a = -1.20 and a = -1.75, which corresponds to the large range of values 
found in the literature due to the lack of constraint on the faint end slope of the LF. 



rest-frame K-band sample, and we conclude that at z < 1 the 
stellar mass density increases by a factor of ~ 2.3. Since the 
global rest-frame ultraviolet emissivity continues to decline, 
merger events should produce little star formation, either via 
minor mergers with e.g. satellite galaxy accretion, or via ma- 
jors mergers between cold gas-depleted galaxies so there is not 
enough gas to efficiently produce new stars. 



At z ~ 1.1 our luminosity densities exhibit a transition in 
the evolutionary trend. In particular, from z = 2 to z = 1.2, the 
FUV-, B-, V- and /?-band X increase, the U- and /-band £ flat- 
ten; and then below z = 1 .2, they all decrease. Only the NUV £. 
follows a continuous evolution since z = 2. Nevertheless, error 
bars at z > 1 .2 are large, and thus the observed transition is still 
uncertain. 



5. Comparison with other surveys 

5.1. Comparison with the CFRS at Ub = 22.5 

The CFRS is an /-selected survey like the VVDS. The CFRS 
galaxy sample consists of 730 /-band selected galaxies at 
17.5 < Iab < 22.5, of which 591 (i.e., more t han 80 percent) 



have re liable redshifts in the range < z < 1 .3 dLe Fevre et al 



19951) . iLiUv etaP (Il996l) estimated the comoving luminosity 
densities of the Universe from the CFRS sample in the rest- 
frame NUV-2800, B-4400, and NIR-IOOOO passbands over the 
redshift range < z < 1 with the cosmology (Q.m, ^a, h) = (1, 
0, 0.5). Here, we compare with the B-band results which re- 
quire very little extrapolation from models, contrary to the 
NUV or NIR data. From flie best estimate of the B-4400 band 



LF (see iLilly et al.L Il995h . the rest-frame B-band emissivities 



directly-observed (correspond to our minimal estimates) and 



estimated were derived. The CFRS estimated X4 



was de- 



rived fitting with a Schechter function the data given by the 



Vmax LF estimator The VVDS estimated X, 



'4400 A 



is derived 
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Fig. 2. The estimated non dust-corrected comoving luminosity 
densities in the rest-frame B-4400 passband as a function of 
log{\ + z) and represented with filled circles. The crosses are 
the minimal comoving densities, derived from the observed lu- 
minosities without extrapolation neither at the bright nor at the 
faint part of the LF. Data are displayed at the mean redshift 
of the galaxies within each redshift bin. Error bars are derived 
from the 2.3cr error contours of the Schechter parameters, ex- 
cept at z < 0.2 where M* is fixed and at z > 1.2 where a is 
fixed. The dot at z = 0.1 is the SDSS local point. The VVDS 
Z?-band luminosity densities increase as (1 H-z)' ''* up to z = 1.1. 
Data are listed in Tab. [1] The open triangles and open circles 
are the DEEP2 and COMBO- 17 data displayed for clarity at 
- 0.05 and -h0.05 in re dshift respectively and taken from Table 2 
in Faber et al. (2005), the open squares are the HDF data from 
Poll et aL (.20031). ari d the open stars are the FDF data from 
Gabasch et al. [ (5oo4 . 



from the integration of the LF estimate with a, M*, and 0* 
as determined with the STY method. The two fits give simi- 
lar results as long as each type of galaxies is visible within the 
redshift range studied (see lllbert et al. , 2004 ). 



The upper panel of Fig. |6] compares the CFRS -C4400 A ^'^ 
the VVDS -C4400 A ^i'^h a magnitude cut at Iab - 22.5 over the 
same redshift range, 0.2 < z < 1.0, and with the same cosmo- 
logy adopted in the CFRS analysis. The CFRS and the VVDS- 
[17.5 - 22.5] X4400 A ^^'"y ^^11 consistent with each other 
(see Tab.|2|. Even though the CFRS contains ~ 4.5 times fewer 
galaxies than the VVDS at the same depth, the CFRS eiTor bars 
are smaller than the VVDS eiTor bars. This is due to different 
procedures used to estimate the uncertainties in the two sur- 
veys. The CFRS ' ad-hoc' uncertainty procedure is described in 
Lilly et al. ( 1996h . and ours are derived from the a - M* error 
contour of 68 percent confidence level of our STY estimate. In 
both surveys, the minimal and the estimated B-band emissivi- 
ties become more discrepant as the redshift increases. This is 



Fig. 3. Comoving non dust-corrected luminosity densities in 
the rest-frame V-5500, R-6500, and 1-7900 passbands. The 
filled circles, crosses, and the dot are the same as in Fig. |2] 
The VVDS V-, R-, and /-band densities increase as (1 -1- z)"^^ , 
(1 + z)" "*^, and (1 -I- z)" "^" up to z = 1.1 respectively. 



Table 2. Comparison of comoving luminosity densities be- 
tween /-selected surveys, the CFRS and VVDS, with the cos- 
mology (Qm, ^a, h)^(l,0, 0.5). 



log £4. 



in WHz-' Mpc- 



Redshift 


CFRS 


VVDS 


VVDS 


range 


[17.5-22.5] 


[17.5-22.5] 


[17.5-24.0] 


Minimal 


z<0.2 






19.50 ±0.03" 


0.20-0.50 


19.61 ±0.07 


19.68 ± 0.01 


19.67 ±0.01 


0.50-0.75 


19.77 ± 0.07 


19.80 ± 0.02 


19.86 ±0.01 


0.75-1.00 


19.74 ± 0.07 


19.71 ±0.01 


19.85 ±0.01 


Estimated 


z < 0.2 


19.30 + 0.10" 




in cq+o.lUn 


0.20-0.50 
0.50-0.75 
0.75-1.00 


19.63 ±0.07 
19.86 ±0.08 
20.05 ±0.13 


19.68!!;- 

i9.88!!;:«i 

20.06^]^ 


19 70+" ''* 
19 90+0 02 

1 Q Q^+0.02 



" With M* fixed to the SDSS ° 'g-LF estimate from Isianton et alJ 
( l2003h (fi-band ~ " 'g-band with less t han 0.05 mag d i fferen ce). 
* Based on SAPM bj-LF estimate from lLovedav et all ( Il992h . 



due to the combination of brighter limiting luminosities sam- 
pled at higher redshifts, and to a lesser extent, a steeper faint- 
end slope of the rest-frame B-band LF The CFRS 0.2 < z < 1 .0 
values are thus found to be very reliable up to Iab - 22.5. 
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Fig. 4. Comoving non dust-corrected luminosity densities in 
the rest-frame U-3600, NUV-2800, and FUV-1500 passbands. 
The filled circles, crosses, open stars, and the dot come from 
the same data as in Fig. |2] The VVDS U-, NUV-, and FUV- 
band densities increase as (1 -i- z)' ''^ , (1 + z)''^**, and (1 + z)^"^ 
up to z = 1.1 respectively. The asterisk is the FUV-hand lo- 
cal point from the GALEX-2dFGRS survey (< z >= 0.055 



Wvder et al.L l2005h . The plain squares are the FUV-band data, 
displayed for clarity at -h0. 02 in redshift, a nd taken from 



Fig. 5. Comoving non dust-corrected VVDS luminosity densi- 
ties in the rest-frame FUV-1500, NUV-2800, U-3600, B-4400, 
V-5500, R-6500, and 1-7900 passbands from bottom to top re- 
spectively, as displayed in Figs. |2] |3] and |4] The solid line 
connects points where the rest-frame band is observed in the 
optical. 



the GALEX-VVDS survey (ISchiminovich et al.L I2005h . The (-^gg Lilly et al 



The CFRS £- 



evolution was found as (1 + 2)3-9±o.75 
19961) . With the same cosmology (Q.m, 



(Baldrvet al..2005) 



5.2. From the CFRS to the 1.5 mag deeper VVDS 

The lower panel of Fig. |6] compares the VVDS X.44QQ ^ 
Iab - 22.5 with the global VVDS at Iab - 24.0 over the same 
redshift range, 0.2 < z < 1.0, and with the same cosmology as 
the CFRS. In the redshift range, 0.2 < z < 1 .0, X4400 A from the 
global VVDS is consistent within the error bars with the one 
from the VVDS-[17.5 - 22.5] (see Tab. ^. Our VVDS local 
value at z = 0.138 is derived by fixing M* to th e SDSS value. 



small open circles are the NUV-hand points from the SDSS /1) = (i, 0, 0.5) and Iab < 24 the VVDS X2800A evolves as 

(1 + 2)2-4±o i. The rest-frame NUV-ISOO is observed in the op- 
tical in the VVDS at z > 0.2, while in the CFRS it is observed 
at z > 0.5. It adds uncertainties due to the extrapolation of tem- 
plates. Thus the NUV CFRS evolution wa s found too steep as 
it was already seen bv lWilson et al.l (120021) . 

In conclusion we find that the B-band luminosity density 
estimated from the VVDS is in excellent agreement with that 
estimated from the CFRS at 0.2 < z < 0.75 (see Tab. EJ. We 
demonstrate that going deeper in magnitude is superior to as- 
sembling larger samples as far as the galaxy sample is not de- 
pendent on cosmic variance and has a well-defined selection 
function. Futhermore, comparing the estimates of the VVDS at 
22.5 mag to the ones of the VVDS at 24 mag shows that the 
error bars determined in the VVDS are well-defined. Indeed 
the estimate at 24 mag is within the error bars of the one at 
22.5 mag since we account for the uncertainty on the LF slope. 
Finally, we confirm that the uncertainties due to template ex- 
trapolation led to the very steep slope of the NUV CFRS lumi- 
nosity density evolution. 

5.3. Comparison with other deep surveys 

In Fig. |2] we added the B-band luminosity densities from 
DEEP2 and COM BO- 17 surveys (values taken in Table 2 in 
Faber et al.L l2005b . from the FORS Deep Fields (FDF) sur- 



(-19. 30 - 51og/!) mag at z = 0.1 (see Table 2 in lBlanton et al. 



2003L with "-'g-band ^ B-band at less than a 0.05 mag dififer- 
ence level), since we do not span the brightest luminosities at 
z < 0.2. In this redshift range, we observe fainter luminosities 
than the SDSS, and we find a s teeper faint-end sl ope than the 



SDSS (see discussion in §5.1 in Ilbertetal., 2005) 



The best fit power law for the estimated -C4400 A of '^^e 
VVDS data gives (1 + ^)2 05±Q Q6 to be co r npare d with the one 
of the CFRS, (1 -hz)2-^2±o.5 (see lLillv et ^111996 ) in cosmology 
(Qm, ^a, h) = (1, 0, 0.5). The steeper slope of the CFRS data 
is due to the combination of the adopted low local reference at 
z = and the high normalization of the CFRS LF at z = 0.85 
due to a poor constraint of a and M* with a sample limited at 
Iab = 22.5. 
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ve y where Schechter pa rameters are taken from Table A.5 



12, e 



7.0 



5.6 



4.6 Gyr 



Gabasch et alj ( 120041) . and from the HDF sur vey where 



Schech ter parameters are taken from Table 2 in iPoli et al 
( l2003b . These surveys at 0.2 < z < 1 are in excelle nt agree- 
ment with the VVDS, except for the IPoU et all (l2003b at 0.4 < 
Z < 0.7 which has a steeper faint-end slope (Aa ~ 0.15) than 
the VVDS, and thus presents a B-band £, higher by 0.1 dex 
than the other surveys in this redshift range. At 1 < z < 1 .2 the 
VVDS is higher by 0. 1 dex than DEEP2 because the faint-end 
slope of the B-band LF is better constrained in the VVDS, a 
deeper sample by ~ 1 mag than DEEP2. We note that both the 
/-selected VVDS and the deeper by 2.8 mag /-selected FDF 
exhibit a drop of the B-band luminosity density at z ^ 1.1. 

Our 1500 A X results are slightly higher than the ones from 
the FDF by 0.1 dex at 0.6 < z < 2. It is likely due to their 
lower fixed {a - -1.07) value than our fixed {a - -1.6) value. 
However, the difference is usually within the error bars. Our 
2800 A £, results are in good agreement with the ones from the 
FDF at 0.4 < z < 1 .2, and in excelle nt agreement with the ones 
from the SDSS dBaldrv et al.[l2005h at z < 0.3. 



6. Luminosity densities per galaxy type at z < 2 

6.1. The VVDS galaxy types 

The large VVDS sample enables us to study the luminosity 
density evolution for different galaxy types. To associate a 
spectral type to our galaxies with a known spectroscopic red- 
shift, we used the best fitting type between UB VRI photometric 
data and a set of SEDs that were lineary interpola ted between 



the four observed spectra of Coleman et al. ( 1980l) , i.e. E/SO 



Sbc, S cd, and Irr, and two st arburst models from the GISSEL li- 
brary (!BmzualandCh^lot,[l993). Then we divided the galaxy 
population into four rest-frame color classes, the elliptical-like 
(type-1), the early spiral-like (type-2), the late spiral-like {type- 
3) , and the irregular- like (type-4) types. We describe in detail 



Zucca et al. ( 2006h the fitting process and the robustness of 



the classification. 

We note that several previous deep surveys were limited 
to the stu dy of two popula tion sub-samples. For instance in 
the CFRS. Imiv et alj (119951) studied blue and red populations, 
simply dividing the galaxy population into two equal num- 
ber sub-samples and corresponding to a Sbc color separation. 
This allowed us to identify little evolution of the red popula- 
tion while the blue population evolves strongly. In more re- 
cent deep surveys, the bimodal rest-frame color distribution ob- 
served at least up to z = 1.5 is being used to defin e an empirical 



separation between red and blue galaxies (K20: iFontana et al 



(2004), DEEP2 andCOMBO-17: iFaber et al.l (I2005h. VVDS 



Ipranzetti et al] (2007)). While the bimodality is clearly ob- 
served for the bright galaxies, the faint population does not ex- 



hibit two contrasting modes (see iFranzetti et al.L 120070 . Also 



the bimodality hides a strong differential, color evolution of 
the bright populat ion (i.e. L > L*) as show n in Fig. 2 of 
Zucca et al. [ (l2006h . Or, the extensive analysis of iDe LapparentI 
( 2003b shows that estimations of LFs based on two color sub- 
samples lacks the necessary discriminatory power for detect- 
ing the variations in luminosity as a function of type which 
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Fig. 6. Rest-frame B-band comoving non dust-corrected lumi- 
nosity densities with the (Om, ^a, h) - (1,0, 0.5) cosmology 
at < z < 1 of the /-selected CFRS and VVDS surveys. 
Values are listed in Tab. |2] For clarity the z > 0.2 data are rep- 
resented by symbols slightly displaced horizontally from the 
center of the redshift bins delimited with the vertical dotted 
lines. The plain square points at z = 0.05, 0.07, a nd 0.1 are 
local points from the SAPM dLovedav et al.lll992h. the 2dF- 



GRS (Norberg et al.L 120021) . and SDSS (IBlanton et al.L 120031) 



respectively. The dashed line shows the best-fit power law and 
its associated uncertainties for the ' LF-estimated' XAdm A of 
the CFRS 17.5-22.5 data derived in iLiUv et al.1 d 19961) . (1 + 
^-)2.72±o.5^ and the sohd line is the one for the VVDS 17.5-24.0 

data,(l +2)2.05±0.06_ 



are traced by the intrinsic LFs. From a t heore tical study us- 
ing semi-analytical models. iKavrraj et al.l (120061) show that the 
red-sequence traces the progenitor set of early-type galaxies in 
terms of numbers and masses for the bright galaxies, but breaks 
down severely at faint (L < L*) luminosities. And thus, using 
the bimodality does not seem to be robust in classifying galax- 
ies in types unambiguously related to passively evolving galax- 
ies on one side and star-forming galaxies on the other Whereas, 
in fact, our selection relies on the complete SED available from 
the multi-band imaging. We note that Seyfert 2 galaxies are in- 
cluded in our sample. We analyse these galaxies with the true 
starburst population in another paper. Nevertheless we do not 
expect a large contribution of the Seyfert 2 population at z < 1 . 

We emphasize that our four types are nicely correlated with 
colors and asymme try-concentration p arameters from HST im- 
ages (see Fig. 2 in lllbert et al.l l2006al) and with spectroscopic 
features (e.g, emis sion line strength, 4000 A break, see Fig. 1 in 
Zucca et al.l 120061) . The four individual LFs were derived bet- 
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Gyr 



■ Elliptiical-liki 
□ Sab-jlike 
A Scd-iike 
O Irr— like 




Redshift 

Fig. 7. Comoving non dust-corrected luminosity densities in 
the rest-frame Z?-passband from early to late galaxy types (see 
details in Section 6. 1). For clarity, data are represented by sym- 
bols slightly displaced horizontally from the center of the red- 
shift bins. Error bars are at Icr. The irregular-like type emissiv- 
ity decreases markedly by a factor 4, while the elliptical-like 
type increases by a factor 1.7 from z - 1.1 to z = 0.2. We 
note that at z < 0.4 the total emissivity is dominated by the 
early-type population. 



ween 0.2 < z < 1 .5 where our sample is essentially complete 
for every type, and Fig.|2]displays the four luminosity densities 
in the rest-frame Z?-band where there is very little extrapolation. 



6.2. Evolution of £ per type 

The evolution of the B-band LP per type and the evolution of 
the fraction of bright (Mfl,„ -5/og(/i) < -20) galaxies has been 



tne traction or brignt (Mfl,„ -3/og(«) < -zu) galaxies nas oeen 
described in detail in Zucca et alT ( 20061) . and consequently, the 
B-band luminosity densities evolve as follows. 

Type-4: irregular, starburst, very blue galaxies. The frac- 
tion of bright galaxies decreases from ~ 35 to ~ 5 percent from 
z— 1 .5 to z = 0.2. The LF undergoes a strong evolution in den- 
sity (negative) and luminosity (negative). There is a strong de- 
crease in volume density by a factor ~ 2 coming from both the 
bright and faint parts of the LF. And thus £, decreases markedly 
from z- 1 .5 to z = 0.2, by a factor ~ 3.5. 

The three other types have an LF which corresponds to a 
mild evolution in density (positive) and luminosity (negative). 

Type-3: late-spirals, star fanning, blue galaxies. The frac- 
tion of bright galaxies decreases from ~ 40 to ~ 10 percent 
from z = 1.5 to z = 0.2, and meanwhile the £ decreases by 
a factor ~ 1.7. To keep decreasing (p* x L* this population is 
strongly faintening by ~ 1.5 mag. X decreases markedly at 



z < 0.4 by a factor ~ 1.5 due to the faintening by ~ 1 mag 
of the LF 

The last two types present an increasing fraction of bright 
galaxies. 

Type-2: early-spirals, post-starburst, red galaxies. The 
fraction of bright galaxies increases from ~ 20 to ~ 35 per- 
cent from z = 1.5 to z = 0.2, and globally X increases by a 
factor ~ 1.6. To keep increasing <p* x L* luminous red galaxies 
must appear at low redshifts since this population is modestly 
faintening by ~ 0.6 mag. Nevertheless there is a significant 
variation at low z. Indeed, £. increases from z - 1.5 to z = 0.6, 
by a factor ~ 1.4, then decreases from z - 0.6 to z = 0.4, by 
a factor ~ 1.3, and finally increases from z - 0.4 to z = 0.2, 
by a factor ~ 1.7. This transition implies that a small fraction 
of luminous type-2 galaxies disappears (~13 percent) while the 
luminosity decreases by a small factor of -0.15 mag. 

Type-1: elliptical, red galaxies. The fraction of bright 
galaxies increases from ~ 0.05 to ~ 55 percent from z = 1.5 
to z = 0.2, and meanwhile X increases continuously , by a fac- 
tor ~ 1.7. That is luminous red galaxies must appear at low 
redshifts to keep increasing (p* x L* since this population is 
faintening by ~0.3 mag only. 

In addition, Fig.|7]shows that the type-1 and type-2 red po- 
pulations dominate the total light at z < 0.4, while at 0.4 < z < 
1 .2 the late-type does. Since we know that red spher oids are the 



majorit y of our type-1 population from our work in lllbert et al 



(l2006ah . we conclude that a dust deficient population is domi- 
nant at z < 0.4. The luminosity density increasing of the type-1 
and type-2 red populations suggests a contribution from merg- 
ing phenomena. Indeed in a downsizing scenario where lumi- 
nous red galaxies are already in place at high redshifts (z » 1), 
and low luminosity red galaxies appear at low redshifts, a flat 
luminosity density would be expected as a function of redshift, 
whilst adding merging would increase the luminosity density 
of the luminous red population as redshift decreases. 

We note also that from z - 0.7 to z = 0.5, X for type-2 
decreases by a factor ~ 1.3, while X for type-1 increases by a 
factor ~ 1.15; this might suggest that <15 percent of ellipticals 
are formed from early-spiral major mergers within this 1.3 Gyr 
period. It might be evidence for mergers between gas-deficient 
bright galaxies. Futhermore, the fact that the type-2 population 
luminosity density increases again from z - 0.5 to z = 0.3 by a 
factor ~ 1 .7 might suggest a density growth due of an evolution 
of type-4 galaxies towards type-2 galaxies. 

The X decreasing of the type-4 population is markedly dif- 
ferent from that of the remaining population. This pop ulation 
is dominated by dwarf galaxies (see IZucca et all l2006l) . The 
evolution of this population supports a downsizing scenario 
where most star formation is shifting to low-mass galaxies at 
z < 1 .2, while the global luminosity density is dominated by 
other galaxy types. 

7. Global FUV luminosity densities at 2.7 < z < 5 

7.1. The VVDS high-z population 

In the high-redshift range, 2.7 < z < 5.0, the rest-frame 1500 A 
corresponds to the observed frames R and / passbands. Thus 
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<z>° 


Redshift range 




fixed 




M;/1500) mag 


^' ilQ-'^ Mpc-^) 


Unweighted 


3.04 


[2.70-3.40] 


-1.6;- 


1.2 ; - 


1.75 


-21.68 ; -21.46 ; -21.76 


6.27 ; 8.02 ; 5.66 


3.60 


[3.40-3.90] 


-1.6;- 


1.2 ; - 


1.75 


-22.52 ; -22.38 ; -22.84 


2.64; 3.45; 1.56 


4.26 


[3.90-5.00] 


-1.6;- 


1.2 ; - 


1.75 


-22.72 ; -22.44 ; -22.86 


0.63 ; 0.92 ; 0.50 


Incompletness-corrected'" 


3.46 


[3.00-4.00] 


-1.4 ; - 


1.2 ; - 


1.75 


-21.38 ; -21.29 ; -21.57 


1.23 ; 1.32 ; 1.00 



"Mean redshift of the galaxy redshifts in the quoted redshift range. 

''We fixed or to - 1 . 6, -1. 2, and -1.75 and the corresponding STY values for M* and 0* are given, respectively. 
^See lPaltani et alj ( l2007h . With a = -1.6. we have M* = -21.49and<^* = 1.11. 



uncertainties on absolute magnitudes related to ^[-corrections 
and galaxy types are small. Our data set consists of 161 red- 
shifts with VVDS quality flags 2, 3, and 4, and 237 red- 
shifts with VVDS quaUty flags 1 in the VVDS-0226-04 field 
(~0.5 deg^). The latter fraction is not negligible since the dif- 
ficulty in determining a redshift is increased at the faintest ap- 
parent magnitudes. The confidence levels for the single high-z 
population are > 50 percent for the flags 2, 3, a nd 4 dataset, and 



~ 45 p ercent for the flags 1, 2, 3, and 4 dataset (ILe Fevre et al 
2005bl) . 



Here, we correct our sample for the target sampling rate 
only (see Section 2). That is, there is no assumption for the 
redshift distribution of the sources that were not spectroscopi- 
cally observed or for which no redshift could be reliably iden- 
tified from the spectrum obtaine d. This is a more restricted ap- 
proach than what was adopted in lllbert et al . (2005), where the 
photometric redshifts were used to obtain a spectroscopic suc- 
cess rate as a function of redshift, and thus to correct further 
for incompleteness. We could not use the same approach here 
because our pho tometric redshifts h ave been thoroughly tested 
only up to z ~ 2 dllbert et al.Ll2006cl) . The correction for incom- 
pleteness has for effect to steepen the slope of the unweighted 
LF, and since a is correlated wi th M*, it produces a brightening 
(see Fig. 4 in lllbert et aill2005b . 



We consider the two following high-z datasets. The stan- 
dard dataset is composed of quality flags 2, 3, and 4 red- 
shifts and the extended dataset is composed of quality flags 
1, 2, 3, and 4 redshifts. It is likely that the true luminosity 
function/density lies between these two cases, assuming that 
the 7 percent of flags over the whole < z < 5 sam- 
ple makes little difference. Using one field of view is likely 
to induce uncertainty due to cosmic variance. Nevertheless 
Sawicki and Thompson ( 2006t) observed an effective area of 
169 arcmin^, i.e. one third of ours, 410 arcmin^, split into five 
fields (the Keck Deep Fields; KDF) and found a quite moderate 
cosmic-variance effect. 



Inl Paltani et al. (2007) we do a comprehensive study about 
the impact of the large uncertain redshift population at 3 < 
z < 4, in using a somewhat different analysis. Here we choose 
to use the same analysis as at z < 2. We note that results are 
qualitatively similar. 



7.2. Estimation of tlie high-z FUV VVDS luminosity 
functions and densities 



For the rest-frame FUV analysis, we fixed a to -1.6. As shown 
in Fig.[T]the uncertainty in a slope a < -1.3 may lead to large 
discrepancies in the estimation of For consistency, we have 
derived error bars of the rest-frame FUV X. from the STY esti- 
mations assuming a = -1.20 and a = -1.75 for the extended 
high-z data set as we did at z < 2. The LF parameters are given 
in Tab.[3] Fig.|8]displays the LF estimations of the standard and 
extended high-z datasets, and the incomplete ness-corrected LF 
estimated at 3 < z < 4 in lPaltani etai](l2007b . We do a detailed 
co mparison for differe nt fixed a values found in the literature 
in ' 



Paltani et al.1 ilOOlV The rest-frame FUV X at 2.7 < z < 5 
are displayed in Fig. |9] together with those at z < 2, and they 
are listed in Tab. |4] 



7.3. Comparison with other high-z surveys 



Fig.|9]displays other rest-frame FUV luminosity densities from 
the literature. At 2 < z < 5, our /-selected rest-frame FUV 
X. are broadly in agreement with the other estimations. At 
z ~ 3, rest-frame FUV X estim ations are found with in 0.1 dex 
(19.47 for the VVDS , 19.58 for Steidel et al. (1999'), 19.43 for 
Gabasch et all (l2004. 19 .5 2 .Arnouts et aL (.2005). and 19.51 
for lSawicki and Thompson ( 20061) ). We note that the lumino- 
sity functions exhibit noticeable discrepancies, in particular at 
the bright end, and furthermore that the faint end slope is no t 



constrained (see detailed comparisons in iPaltani et al.L l2007n . 



As at z < 2, our high-z luminosity densities are slightly higher 
than those estimated from the FDF (due to their lower adopted 
slope of a - -1.07) except at 3.9 < z < 5. The decrease from 
z = 3.4 to z = 5 is smaller in the FDF (0.25 dex) than in the 
VVDS (0.5 dex), however our error bars are compatible with 
the FDF result. This may suggest that our unsuccessful spec- 
troscopic identification rate has a strong effect in our highest 
VVDS redshift bin, while our data points at 2.7 < z < 3.9 
appear fully consistent with the FDF. 
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8. The FUV luminosity densities at < z < 5 

8.1. The global shape 

Fig.|9]displays the rest-frame FUV luminosity densities within 
the redshift range < z < 5. We find that the global luminosity 
density increases by a factor ~ 3.5 from z - 0.05 to z = 1.2, 
by a factor ~ 1.2 from z = 1.2 to z = 3.4, and decreases by a 
factor ~ 0.3 from z = 3.4 to z = 5. Futhermore the evolution at 
1.1 < z < 3.0 might be more complex than a modest increase. 
Indeed, even though our large eiTor bars cannot exclude an in- 
crease or a plateau, our data points taken at face value exhibit 
a decrease by a factor ~ 1.3 from z - 1 to z = 2, and an in- 
crease by a factor ~ 1.4 from z - 2.0 to z = 3.4. One could 
mention dust attenuation effects since at these redshifts, targets 
are selected from their rest-frame UV, and thus a fixed /-band 
flux cut could miss a non-negligible fraction of highly dust en- 
shrouded targets. However, UV luminosity densities at z ~ 3 
are larger than at 1 < z < 2 while they should be even more 
dust affected since selected at shorter ultraviolet wavelengths. 

The 1500 A luminosity density as a function of redshift 
from the FDF (iGabasch et al.ll2004l) presents a similar shape to 
that derived from the VVDS. We stress that both datasets, the 
FDF and the VVDS, span the redshift range 0.5 < z < 5 within 
one single survey and have the same unique /-band selection 
criterion. We note that the FDF used photometric redshift tech- 
niques with NIR photometry, which gives reliable redshifts at 
1 < z < 2. The VVDS and the FDF are complementary; the 
VVDS is ~ 50 times larger in surface than the FDF and it con- 
sists of spectroscopic redshifts, while the FDF goes 2.8 mag 
deeper than the VVDS using photometric redshfits. Both esti- 
mates are in agreement, in particular at 1 < z < 2. 

We note tha t the empirical models from 
Perez-Gonzalez et al.l (l2005h also show a decrease in the 




STY (a fixed 



-20 -1 
Mj3(1500A) 



5 log h 



Table 4. Comoving FUV- 1500 luminosity densities at 2.7 < 
z < 5 of the VVDS for the extended high-z dataset. 



Fig. 8. Unweighted non dust-coiTected luminosity functions in 
the rest-frame FUV- 1500 passband estimated with a fixed to 
-1.60. The circles and the solid line represent the V„,a j and the 
STY estimates for the extended dataset where the number of 
galaxies is given within parenthesis next to the redshift range. 
The vertical dot-dashed line corresponds to the faint absolute 
magnitude limits considered in the STY estimate (see § 3.2). 
The hashed area is limited by the LFs derived with a fixed to 
-1.20 and -1.75. The dotted lines represent the STY estimates 
with a fixed to - 1 .60 for the standard dataset. The standard 
dataset LF is below the extended dataset LF since it does not 
include quality flags 1 (see §7.1). In the top-left panel, the long- 
dashed line is the incompleten ess-corrected dataset FUV- 1500 
LF at 3 < z < 4 derived as in IPaltani et al.l (l2007h . with fixed 
a= -1.4. 



"Mean redshift of the galaxy redshifts in the quoted redshift range. 
''Elapsed time in Gyr in the redshift bin. 

■^The STY LF estimators are derived using a fixed a = -1.6 and error 
bars are derived from the STY estimations assuming a = -1.20 and 
a = -1.75 as done at z < 2 (see Tab.[Tll. 

''in IPaltani et aP l l2007h . we use a = -1.4. Here, we use the same 
minimal and maximum values as at z < 2 for a, i.e. -1.2 and -1.75. 



<z>" 


Redshift range 


St" log £1,00 A 

WHz ' Mpc-3 


Unweighted minimal 


3.04 


[2.70-3.40] 


0.6 18.82 


3.60 


[3.40-3.90] 


0.3 18.95 


4.26 


[3.90-5.00] 


0.3 18.89 


Unweighted estimated 


3.04 


[2.70-3.40]'' 


0.6 19.47!^:;',' 


3.60 


[3.40-3.90]'- 


0.3 ^9AT°,fg 


4.26 


[3.90-5.00]" 


0.3 18.89+°;^ 


Incompleteness-corrected estimated 


3.46 


[3.00-4.00]'' 


0.5 19.47!!J;f,' 



cosmic SFR density between z = 1.4 and z - 2.2, varying 
from 10 to 29 percent depending on the model but with 
luminosity evolution solely. It is interesting to see that their 
model including a combined luminosity plus number density 
evolution does not exhibit a decrease. Rather than excluding a 
number density evolution it indicates that the global emissivity 
is dominated by the luminosity evolution. We go into more 
detail about this effect in the next section. 

We emphasize that our /-selected sample misses the very 
red (I-K) > 2.6 galaxy population at 1 .2 < z < 2.5 . This popu- 
lation concerns the most massive star forming dusty galaxies. 



With a 



-l.60^°jlwe find log £: 



19 64+" 20 



as detailed in our /iT-selected sample described in lPozzetti et al 
(12007,1 . 



8.2. The intrinsic luminosity dependency 

The luminosity density is dominated by the evolution of the 
bright population. Thus we have integrated the LF from the 
three following absolute magnitude limits: Mab(1500 A) < 
-19 mag, Mab(1500 A) < -20 mag, and M^b(1500 A) < -21 
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Fig. 9. Comoving non dust-corrected FUV luminosity densi- 
ties from z = to z = 5. The plain circles represent the 
1500 A VVDS data (0.05 < z < 5) (see values and error 
bars in Tabs. [T]&|4]|, the plain triangle r epresents the 1500 A 



VVDS data at 3 < z < 4 derived as in iPaltani et al.1 (l2007h . 



and the plain squares represent the 1500 A GALEX-VVDS 
data at z < 1, plo t ted at -1-0.02 in redshift for clarity, from 
Schiminovich et al.l(l2005b. Othe r data are the 1500 A GALEX- 
2dFGRS from Wvder et al. (2005) (asterisk), the 1500 A HDF 
data from Arnouts et al.l (120 05) (crossed circles), the 1700 A 
KDF data from ISawicki and Thompson (20060 (open rhom- 
bus), the 1500 A FDF data from iGabasch et atl (|2004 (open 
stars), and the 1700 A data from Iwata et alT i 2007 ) (crossed 
square). 



mag. According to our LFs, these magnitude cutoffs corre- 
spond about to L > L* galaxies at 0.1 < z < 0.6, 0.6 < z < 2, 
and 2 <z < 3.4 re spectively. At z ~ 0.05, M* (1500 A) 



-18 

mag dWvder et al.L [2005). We choose a fixed cutoff rather than 
a luminosity evolving cutoff to make a comparison over a long 
time baseline because this does not depend on the modeling 
of the luminosity evolution of the LF. The change of M* with 
redshift implies that at z - 3 the most luminous galaxies were 
forming stars at a rate ~5, and ~7 times higher than at z - 0. 15, 
and z - 0.05 respectively. We have applied the same cut- 
offs t o the LF integration of the 1500 A HDF ( Arnouts et"aLl 
2005'). the 1700 A KDF ("Sawicki and Thompson', '2006^ and 
the 1700 A data point of Steidel et al. (1999) . Fig.fTOldisplavs 
these different luminosity densities for the bright population. 
We also produced the rest-frame FUV-band luminosity densi- 
ties for three range of luminosities, Mab(^5Q0 A) > -19 mag 
(dwarf and intermediate population, < 0.1 L* „ i r- i ), 

^ i- i- ' Lyman Break Oalaxy^' 

-21 < Mab(J500 A) < -19 (luminous population), and 
M/iB(1500 A) < -21 (very luminous population, > L'g^), as 
shown in Fig.fTTl 



8.2.1. Detailed results 

We observe that on average the rest-frame FUV-band lumino- 
sity density of the (2.7 < z < 5) population is significantly 
brighter by a factor at least ~ 6 than the (0.2 < z < 1 .4) popula- 
tion. This implies a transition phase for the population dominat- 
ing the FUV emissivity. Another transition phase is observed 
at z < 0.2 with a steep decline by a factor of at least ~ 10 
of the bright population luminosity density. Concentrating on 
the most luminous galaxies, and taking the values of X at their 
face value, we can tentatively identify five phases from z = 5 
to z = as follows. 

(a) From z = 5 fo 3 .4, the emissivity due to galaxies brighter 
than Mab(1500 A) < -20 and -21 mag increases by a factor 
of ~ 3 and 2.5 respectively. The whole population, from the 
dwarf to the very luminous galaxies, sees its FUV luminosity 
increasing. This corresponds to a very active phase in terms of 
newly formed stars in every galaxy. 

(b) From z = 3.4 fo 1.4, the emissivity due to galaxies 
brighter than M^b(1500 A) < -20 and -21 mag steadily de- 
creases by a factor ~ 6 and 25 respectively. This ~ 3Gyr phase 
corresponds to the progressive, but relatively quick, drop of star 
formation in the most luminous galaxies which are formed ear- 
lier than z = 3.4, and whose contribution to the global FUV 
lumin osity densit y beco mes less important toward lower red- 
shifts. Lotz et al. I (|2006') observe that the distribution of the 
HST rest-frame FUV morphologies at z ~ 1.5 are similar to 
the ones at z ~ 4 with an identical fraction of major-merger 
candidates. This suggests that this drop is due to the same lumi- 
nous galaxy population, and that it is not caused by a decrease 
in the number of mergers. Thus from z - 3.4, the most FUV 
luminous galaxy population has suddently finished its stellar 
mass assembly and active star forming phase, while the less 
FUV luminous population is still very active and becomes the 
dominant star forming galaxy population until nowdays. 

(c) From z = I A to 0.6, the emissivity due to galaxies 
brighter than Mab{^50Q A) < -20 and -21 mag increases by 
a factor ~ 2 and 4 respectively. This corresponds to Fig. 2 in 
Zucca et al. ( 20061) . where the fraction of the bright population 
(i.e. L > L*) of early-type (E- and Sab-like) galaxies increases 
while the one for late-type (Scd-, Irr-like) galaxies decreases 
even though the latter still dominate in terms of numbers and 
also in emissivity (see Fig.|2ll up to z ~ 0.7. This phase could 
correspond to merger events; indeed very disrupted HST mor- 
phologies in the very luminous galaxy population are generally 
observed. In the CDFS, we observe mergers with i n the e arly- 
type population as shown in Fig. 2 of lllbert et al. (2006a) and 
also that the volume density of red, bright bulge-dominated 
galaxies increases by a factor ~ 2.7 from z ~ 1 to z ~ 0.7 . 
Furthermore from z = 1.5 to z = 0.6 we find in lllbert et al 



(.2006b ) that rest-frame B-band X for galaxies lying in over- 
dense environments (at a scale of 5/1"' Mpc) increases by a 
factor 1 .2, while the one for those in under-dense environments 
continuously decreases by a factor of 2.3. Since the most lu- 
minous galaxies are usually found in dense environments, our 
results favor the merging event to build-up early-type galaxies. 
The exact percentage of merger events is still hotly debated. 
We note also that it corresponds also to the phase where the 
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less luminous population has reached its peak of star forma- 
tion activity at z ~ 1.1. Another point is that at 1.5 < z < 2.5, 



our / -selection is missing the reddest galaxies dPozzetti et al 



20071) . the dusty a nd mas s ive so urces as discovered in infrared 
surveys (i.e., e.g. iDaddil (|2004 ). and undergoing strong star 
forming and dusty phase. This population would be an excel- 
lent contributor to populate the FUV luminous population at 
z < 1 .4 once the dusty phase linked to intense burst does not 
dominate anymore. 

(d) From z - 0.6 to 0.2, the emissivity due to galaxies 
brighter than Mab{^50Q A) < -20 and -21 mag decreases by 
a factor ~ 1.3 and 1.5 respectively. The decrease of the bright 
part of the emissivity corresponds to a phase entirely domi- 
nated by the early-type galaxies for which the star formation 
decreases passively. 

(e) From z - 0.2 to 0.05, the emissivity due to galaxies 
brighter than Mab(1500 A) < -20 and -21 mag decreases 
by a factor ~ 100 and 25 respectively; the bright part of the 
global FUV luminosity de nsity abruptl y drops. Analysing the 
< z < 0.3 CFRS spectra. iTresse et alj (fl996.) found that only 
54 percent of galaxies exhibit strong star formation with Ha 
and HjS both in emission and 15 percent have Ha and Hy6 both 
in absorption. The remaining 31 percent is the intermediate 
population which shows spectral features from both star burst 
and quiescent stellar evolution. Therefore star formation seems 
to have slowed down or stopped in a significant fraction of 
galaxies in the nearby Universe. This population might corre- 
spond to the descendants of L* galaxies at 0.2 < z < 1 .4 which 
when they fade in luminosity, stop dominating the rest-frame 
FUV-band luminosiy densities. We note that the local sudden 
drop might be combined wit h the local un derdensity seen in op- 
tical re dshift surveys tZucca et al. ( 1997 ). iTresse and Maddoxl 



In resume, we observe a first major crisis for produc- 
ing efficiently new stars, lasting ~ 3 Gyr, from z ~ 4 to 
z ~ 1, which involves the most massive and luminous galaxies 
(Mab(1500 A) < -21, M. > lO^M©). While gradual fading 
of the global population starts at z ~ 1, we observe a second 
major crisis which started ~ 3 Gyr ago when star formation 
was progressively stoping in intermediate mass and luminosity 
galaxies (-20 < Mab(1500 A) < -19, 10'^ < M. < IO^Mq). 
The period from z - 1 .4 to z = 0.6 is highly intricate due to in- 
tertwining populations that follow different evolutionary path. 
Nevertheless the combined effects of decreasing both luminous 
and intermediate population is to make the global luminosity 
density decrease faster at z < 1.2 than at 1.2 < z < 3.9. 



8.2.2. Evolution of L > L* galaxies 

According to the Sc hmidt law of star-forming galaxies 



jKennicutt et al.L Il998l) . the SFR is scaled to the cold gas den- 



sity to the 1 .4 power. Thus, once the gas reservoir is exhausted, 
star formation ceases. The old, most luminous and massive 
galaxies have exhausted their gas reservoir during their early 
intense star formation z > 3.5 (see previous §), and since then 
they undergo passive evolution as star formation ceases. This 



creates excellent dry candidates, that is cold gas-depleted to 
prevent new star formation. 

Creation of new galaxies occurs as the threshold ampli- 
tude for formi ng bri ght galaxies decreases as described in 
Marinoni et^ (l2005b from the VVDS data. That is, the typ- 
ical L* of the population created at a given z will decrease with 
decreasing redshift. This implies that younger, less massive and 
less luminous L* galaxies continue to efficiently form stars with 
a large reservoir of cold gas at 0.2 < z < 3.9. This intermedi- 
ate galaxy population presents a peak of SFR at z ~ 1 . And at 
z < 0.2 these galaxies appear to have also exhausted their gas as 
suggested by the 1500 A X. drop, and start to evolve passively. 
The gas-exhaustion would favor the evolution of morphologies 
toward early-type galaxies. 

In resume, the global FUV luminosity density decreases 
faster at z < 1.2 than at z < 3.9 since it combines the decrease 
from both luminous and intermediate populations. During the 
phase 1.2 < z < 3.9 it undergoes a decline of 0.06 dex, while 
for the phase 0.2 < z < 1.2 it is 0.45 dex. The small de- 
crease corresponds to the phase where the dwarf and interme- 
diate galaxy population increases its SFR, while the luminous 
population start to undergo its SFR decline. In Fig.|9] the global 
luminosity density might decrease by ~ 0.1 dex from z = 3.9 
to z ~ 2, then increase by ~ 1.2 dex from z ~ 2 to z = 1.2. 
Nevertheless, our large error bars do not exclude a smooth de- 
crease by 0.06 dex from z = 3.9 to z = 1.2. 

At z S 4, both processes, i.e. dry mergers toward decreasing 
redshift and morphologies evolving toward early-type galax- 
ies, might contribute to an increase of the bright early-type 
population by a factor ~ 10 to reach ~ 55 percent of the to- 
tal population from z - 1.5 to z = 0.2, while the early-type 
po pulation undergoes a passive luminosity evolution as shown 
by Zucca et aP (2006|). It i s in ag reement with semi-analytical 
studies, e.g. lKaviraj et alJ (120061) find that less than 50 percent 
of the stellar mass which ends up in ear ly-types today is actu- 
ally in early-type progenitors at z ~ 1, or lDe Lucia et al. I (l2006l) 
find that 50 percent of local elliptical accrete half of their stellar 
mass out z ~ 0.8. 

9. The history of the star formation rate density 

We derive the SFR densities hich are not dust corrected using 
the rest-frame 1500 A luminosity densitie s of the VVDS from 
z = to z = 5. The SFR calibration of iMadau etaP (Il998l) 
y ields SFR^ 8 x 102^L(1500 A/erg s"' Hz"!) Mq yr^^ for 
a ISalpeteii ( Il955h 's IMF including stars from 0.1 to 125 so- 
lar mass. We recall that at < z < 1, our 1500 A rest-frame 
band spans ultraviolet wavelengths shorter than 3000 A, that 
is a non-observed wavelength range with the optical bands. 
Nevertheless, we checked that our results are fully consistent 
with the rest-frame 1500 A GALEX-VVDS data (see details in 
Section 4.2). We present our SFR density in Fig.[T2l 

9.1. Dust obscuration 

Our FUV data are not dust corrected. In the following, we com- 
pare them to data which account for dust at z < 2.0 to estimate 
the global correction needed to recover the total SFR density. 
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Fig. 10. Comoving rest-frame FUV non dust-corrected lumi- 
nosity densities from z = to z = 5 for three bright populations 
defined as Mab(1500 A) < -19,-20, and -21 mag. Symbols 
are the same as in Fig. |9] Error bars correspond to the fixed 
slopes, Q- = -1.2 and -1.75 at 0.2 < z < 5, at 0.05 < z < 0.2 
to the fixed M* +0.1; those which are not visible are smaller 
than the size of the data point, i.e. < 0. 1 dex. We do not plot the 
local values for M < -21 which have low values, at z = 0.14 
and 0.055 Log £ = 12.45 + 0.55 and 12.10 Watt Hz"' Mpc-^ 
respectively. In each panel, the solid line connects the VVDS 
points of the global FUV luminosity density as displayed in 
Fig- HI and the vertical lines delimit the approximative redshift 
range in which the bright population is more luminous than the 
characteristic FUV luminosity of the whole population. 

That is, we us e the Ha nebular e mission line dust-corrected 
data taken from Tresse e t al. 1 (12002^ and the 1 2yum m id-infrared 
continuum data from jPer ez-Go nzalez et al.l (l2005h . These re- 
sults are both described with an evolution proportional to ~ 
(1 -I- z)"*, while the evolution of 1500 A X in the VVDS is best 
described by ~ (1 -H z)^ (see Section 4. 1). Assuming an average 
dust coiTection for the whole field galaxy population, the dust 
obscuration at 1500 A is ~ 1.8-2 mag from z = 2 to z = 0.4. 
And from z - 0.4 to z = 0, it becomes smaller from ~ 2 down 
to ~ 0.9 - 1 mag. 

The apparent dust transition at z < 0.5 corresponds exac- 
tly to the change of the dominant population as seen in Fig. [T] 
where the early-type galaxy Z?-band emissivity starts to domi- 
nate below z = 0.4. The early-type population is known to be 
dust deficient, and thus the needed amount of dust obscuration 
is less strong as redshift decreases. Our argument assumes that 
the increasing early-type dominant population in the B-band 
emissivity holds in the f/V-band emissivity. From z = 0.5 to 
z = 2, the B-band emissivity is dominated by the late-type 
population, which form stars from gaseous nebulae, and thus 
is attenuated by a constant factor. From z = 2 to z = 4, the 



Fig. 11. Comoving rest-frame FUV non dust-coiTected VVDS 
luminosity densities from z = to z = 5 for three luminosity- 
class populations defined as MyiB(1500 A) > -19 (dwarf 
and intermediate-luminosity galaxies, < 0.1 Llbq), -21 < 
Mab(1500 A) < -19 (luminous galaxies), and M^b(1500 A) < 
-21 (very luminous galaxies, > L*g(, ). We used and connected 
the same VVDS dataset as displayed in Fig.fTOl 



star-forming high-mass galaxies still dominate the emissivity, 
and thus we expect a constant dust attenuation of 2 mag. We 
note that near-infrared surveys miss the faint blue galaxies at 
z < 1 .2, while they include the extremely red galaxies at z > 1 .2 
as shown in Pozzetti et al. ( 2007). Still the 12//m data from 
Perez-Gonzalez et al. 120051) exhibit a similar shape than the 
FUV data, which reinforces the no evolution of the dust content 
as shown by Bell et al. (2004) at z < 1 and Reddy et al. (2006) 
at z > 1 . The dust correction might change again at z > 4 where 
the by-products of the physical processes to actively form star 
in massive galaxies will dominate through AGN and SN feed- 
back. 

We will make detailed dust obscuration estimates at < 
z < 2 in future papers using the VVDS-GALEX and -SWIRE 
data. 



9.2. Comparison to simulations 



In Fig . [121 display the si mulatio ns from Somerville et al 



(l200lb . Icrotonetal 



(l2006l) . 

Comparing simulations with observational 
might be severely affected by the assumed 
dust obscuration. These three simulations use 
population synthesis m odels of Bruztial and Chariot ( 1993 ). 
Naga mine etaP ilOOdt) use s the IMF of IChabrieil ('2003"). while 
the oth ers use the one of Salpeterl (Il955b . The SFR calibra- 
tion of iNagamine et al. I (l2006l) yields a conversion factor of 



and iNagamine et al 



(l2006l) . 

measurements 
IMF, and the 
the stellar 
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Fig. 12. Star formation rate densities as a function of red- 
shift. Symbols (plain large circles, plain large squares, aster- 
isk, crossed circles, open rhombus, open stars, crossed square) 
and error bars are the same as in Fig. |9] and they represent the 
SFR densities with no reddening correction. The dotted, long- 
dashed line connects the reddeni ng corrected Ha data (small 



plain circles) fro m the CFRS (ITresse and Maddoxl 11998 



Tresse et al.L 120021) . The dotted, short-dashed line connects the 



12/xm data (small plai n squares) from iPerez-Gonzalez et al 



(2005). Models from Somerville et al.1 (1200 lb . iCroton et al 



(2006) and Nagamine et al. ( 20061) are represented with the 



dashed, the dotted, and the short-lo ng dashed lines respectively. 
Everything is for a lSalpetei ( 1955 1) IMF. The derived Apuv dust 
obscuration between the VVDS data not dust corrected and the 
12yum data at z < 2 is represented in the inset. 



1.24 X 10^** for a IChabrieil (l2003h IMF including st ars from 
0.01 to 100 solar mass. It differs by ~0.2 in log from a lSalpeter 
(119551) IMF, and t hus, for consistency, w e add this factor to 
the simulations of Nagamine et al. ( 20061) . The dust prescrip- 
tion differs in a complex manner in each simulation. We do not 
attempt to homogenize the simulations in terms of dust obscu- 
ration. 

The simulations all exhibit a peak of the SFR density at 
z ~ 4- with a smooth decrease up to z = 0. At a quick glance, 
they do not seem to go through the observational points, either 
dust corrected or not. However, at 2 < z < 3.9, the simula- 
tions are close to our FUV data taken at face value if we as- 
sume a constant dust attenuation of ~ 2 mag as discussed in 
the previous paragraph. And, at < z < 0.2, the simulations 
are in good agreement with the dust coiTected data (Ha and 
12 ^m), or with FUV dust -corrected data by 1 mag. We note 
thatlNagamine et al. I (l2006h used the dust extinction factors of 
Steidel et al.1 d 19991) . i.e. 1 mag for z < 2 and 1 .7 mag for z > 2. 
Still, at 0.2 < z < 2, the simulations do not reproduce the 
observational data. At z < 0.2 and z > 2, the global FUV emis- 



sivity is dominated by luminous, massive, large galaxies, while 
at 0.2 < z < 2 it is dominated by the intermediate population. 

10. Conclusions 

We studied the first epoch VVDS data purely /-selected at 
[17.5-24] mag in (AB). The sample is unique in the sense that 
it goes deeper than previous /-selected spectroscopic samples, 
it has a well-defined selection function and it has enough data 
to study sub-samples. Within a single survey, we trace the evo- 
lution of the galaxy population dominating the total light at 
different redshifts all the way from z = 5 to z = 0. The main 
results of our comprehensive study are summarised below. 

i) To study the luminosity density evolution, observing deeper 
in flux is superior to assembling large samples as far as the 
galaxy sample is not dependent on cosmic variance and has 
a well-defined selection function. Nevertheless, the survey 
field must be large enough to identify the rare bright objects 
required to constrain the bright end of the LF. 

ii) The luminosity density evolution is substantially 
wavelength-dependent since the rest-frame passband 
luminosity is related more or less directly to very different 
stellar populations. We find that the non dust-corrected X. 
evolve with time over 0.05 < z < 1.2, as X k (1 h- z)" 
with X = 2.05, 1.94, 1.92, 1.14,0.73,0.42, and 0.30 in the 
FUV-1500, NUV-2800, U-3600, B-4400, V-5500, R-6500, 
and 1-7900 passbands, respectively. Although error bars are 
still large, we note that most luminosity densities exhibit 
a transition at z - 1.1 in the evolutionnary tendency. Over 
the last 8.5 Gyrs, the SFR-related FUV-band luminosity 
density drops by a factor 4. Recent merger events should 
produce little star formation. 

iii) The B-band luminosity density evolution is strongly type- 
dependent. From z - 1.2 to z = 0.2 the irregular-like, 
and Scd-like type galaxies decrease markedly by factors 
3.5, and 1.7, respectively while the elliptical-like, and Sab- 
like types increase by factors 1.7, and 1.6, respectively. The 
late-type galaxy population undergoes a downsizing sce- 
nario where most star formation is shifting to intermediate- 
and low-mass galaxies at z < 1.2. The early-type galaxy 
population suggests a little contribution from merging or 
simple rest-frame color evolution phenomena, but it ev- 
oles passively in terms of luminosity, which suggests a very 
early formation stage (z » 4) for the majority of this red 
population, j 

iv) The SFR density as seen in the rest-frame FUV-1500 lumi- 
nosity density without dust correction undergoes the sev- 
eral following up-and-down phases as redshift decreases, 
and is strongly luminosity dependent. From z = 5 to 
z - 3.4, it increases by at most a factor ~ 3.5, and it cor- 
responds to the end of the mass assembly of the most lu- 
minous galaxies, M/ib(1500 A) < -21 (L>L*gj,), which 
present a peak of SFR at z - 3.5. From z = 3.4 to z = 1.2 
it decreases by a factor 1.2. Notwithstanding the large error 
bars this small decrease might be the result of two oppo- 
site processes; a decrease by a factor ~ 1.4 from z = 3.4 to 
z = 2 due to the fading of the giant galaxy population by a 
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factor 40, and an increase by a factor ~ 1 .3 from z = 1 .2 
to z = 1.0 due to the shift of the star formation activity 
towards the less luminous (-21 < M^b(1500 A) < -19) 
galaxy population which presents a peak of SFR at z ~ 1 . 
From z - 1.2 to z = 0.05 it declines steadily by a factor 4, 
since this phase undergoes the fading of both the giant and 
the intermediate galaxy populations. Nevertheless it hides a 
strong SFR drop by a factor 100 at z < 0.2 of the intermedi- 
ate galaxy population. This phenomena might be combined 
to the local underdensity seen in optical redshift surveys 
(see IZucca et alj (Il997l) . iTresse and Maddoxl (Il998l) ). and 
thus the local Universe might not be an excellent reference 
local point to trace back the evolution of the global galaxy 
population. Our observed global evolution does not seem 
to be in agreement with a continous smooth decrease from 
z ~ 2 to z - as predicted by the simulations, 
v) Comparing our SFR FUV-derived densities with mid- 
infrared or Ha SFR-derived densities, we find that at 0.4 ^ 
z 5 2 the FUV is obscured by a constant factor of ~ 1.8-2 
mag, and at z < 0.5 it is progressively less obscured down 
to ~ 0.9 - 1 mag. In parallel, we find that from z = 0.4 to 
z - 0.05 the Z?-band luminosity density is more and more 
dominated by the early-type (E/Sab-like) galaxy population 
which is known to be dust deficient. Futher analysis com- 
bining VVDS-GALEX-SWIRE data will refine this result. 
We conclude that the old, most luminous and massive 
galaxies have exhausted their cold gas reservoir during their 
early intense star formation which has occured in the early 
Universe at z » 4, and since z - 3.5, i.e. ~ 12 Gyrs, they 
undergo passive evolution as star formation cease. This cre- 
ates excellent dry merger candidates, that is, in which gas 
has been sufficiently depleted that it prevents new star for- 
mation. The small level of contribution of dry-mergers re- 
mains to be quantified. Younger, less massive and less lu- 
minous L* galaxies continue to efficiently form stars with 
a large reservoir of cold gas up to z = 0.2. And at z < 0.2 
these galaxies appear to have also exhausted their gas sup- 
ply as suggested by the FUV-hand luminosity density no- 
ticeable drop, and start to evolve passively, and might be 
combined with a local underdensity. This picture is con- 
sistent with the dow nsizing scenario for the star formation 
jCowie et al. ■ Il996l) . We recall that the luminosity and SFR 
densities are mainly dominated by the L*(p* galaxies rather 
than the dwarf population which is usually undergoing den- 



vi) 



sity evolution (lllbert et all l2005h . At z < 3.5 dry merg- 



ers and morphologies evolving towards early-type galax- 
ies might contribute to increase the number density of the 
bright early-types in n iaintaining a passiv e luminosity evo- 
lution as observed by Zucca et al. (l2006f) at < z < 1.5, 
and also to i ncrease the emissivity in over-dense regions as 
observed by lllbert et aP (l2006b,) . 
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